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^\ 1 Introduction 



Jacobi matrices deserve continuous attention in operator theory, as they play an important role in various areas 
such as numerical analysis, orthogonal polynomials, continued fractions. They also find applications in mathe- 
matical physics: they are discrete analogues of second-order linear differential operators of Schrodinger type on 
the half-line, and appear in integrable systems and the theory of random matrices. 

An important domain of application concerns the theory of orthogonal polynomials of one argument x - 
denoted p„(x) below. Indeed, it is well-known that every sequence {pn}n=o^ satisfies a three-term relation 

XPn = bnPn+l + ClnPn + CnPn-1 , ?^ = 1, 2, ... (1) 

where po = 1, cq = by definition, and an,bn, Cn G C. A Jacobi matrix J being a tridiagonal matrix of either finite 
or infinite dimension, in the infinite dimensional case the non- vanishing coefficients a„, 6„, c„ yield to consider the 
spectral problem associated with the sequence {pn}n=o''- 



Jpn{x) = XPn{x) with 



/ ao ci \ 

feo fll C2 

bi a2 ■ 



\ / 



(2) 



where po{x) , pi{x) , p2{x) , ... form a basis of C[x]. In the finite dimensional case, the main submatrix of J which 
is called the truncated matrix and denoted J*^^^ is usually introduced. If it is of size N + 1 x A^ + 1, then the 
zeros oi pn+i{x) coincide with the A^ -I- 1 distinct eigenvalues of J^^'. In other words, pn+i{x) is proportional to 
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the characteristic polynomial of J^^' and x takes discrete values. As a consequence, the corresponding system of 
orthogonal polynomials is defined on a discrete support. 

Among the families of orthogonal polynomials in one variable, the only ones satisfying ^ and a second-order 
differential equation are known to be the Jacobi, Hermite, Laguerre and Bessel polynomials, as shown in a theorem 
of Bochner [1]. This led the authors of [2] to derive a "g— version" of this theorem using an operator identity of 
independent interest. For this family of orthogonal polynomials satisfying ^ and a second-order g— difference 
equation of the form 

a{y){Pn{s{qy)) - Pn{s{y))) + b{y){pn{s{q^'^y)) - Pn{s{y))) = OnPn{s{y)) for n > , (3) 

where the argument is either x = s{y) = {y + y^^)/2 or x = s(j/) = y, the only solutions to ^ and ([3]) are given 
by the Askey- Wilson polynomials [3J or the big q— Jacobi polynomials, respectively [2]. 

In view of the tridiagonal structure of the l.h.s. of above equations, it is natural to study the algebraic structure 
associated with a pair of Jabobi matrices having a spectral problem of the form ([2]), ([3]). In particular, there 
has been different ways to characterize the Askey- Wilson polynomials and related (7— hypergeometric functions 
from this algebraic point of view. For instance, let us mention 4j where finite dimensional representations of the 
quantum group Uq{sl2) arise, and more generally [5] where finite dimensional representations of the Askey- Wilson 
algebra ^14^(3) are introduced and studied in details. 

In the last few years, a unified algebraic framework based on the concept of Leonard pair has been introduced 
in connection with the orthogonal polynomials of the Askey scheme. Roughly speaking (details can be found in 
[6j), a Leonard pair A, A* acting on a finite dimensional representation V is such that A, A* are diagonalizable 
where the matrix representing A* (resp. A) is irreducible tridiagonal in the basis which diagonalizes A (resp. 
A*). In particular, given a Leonard pair there exists a sequence of scalars /3,j,j*,g,g*,uj,ri,if taken from an 
arbitrary field K such that [7] 

A^A* ~(3AA*A + A*A^-j{AA*+A*A)-gA* ^ j*A^ + ujA + tjI , 
A*^A-pA*AA*+AA*^--/*{A*A+AA*)-g*A = jA*^ + u;A*+7j*I . (4) 

The sequence is uniquely determined by the Leonard pair provided the dimension of V is at least 4, and the 
equations above are called the Askey- Wilson relations (see also |5i). Remarkably, all known examples of Leonard 
pairs for (3 — q + q^^ and q ^ —1 are related with orthogonal polynomials of the Askey scheme [8] in the 
following sense: the entries of the transition matrix (sometimes called the overlap coefficients) relating the two 
"dual" basis which diagonalize A, A*, respectively, can be expressed in terms of one of the following orthogonal 
polynomials ((7— hypergeometric functions): Racah (4-^3), Hahn and dual Hahn (3F2), Krawtchouk (2-F1), g-Racah 
(403), g-Hahn and dual g-Hahn (302), g-Krawtchouk - classical, affine, quantum, dual - (2^1). 

Given the connection between Leonard pairs and the polynomials of the Askey-scheme, several features led to 
consider a more general object called a Tridiagonal (TD) pair 9J. From an algebraic point of view, tridiagonal 
pairs play an important role in the representation theory of the tridiagonal algebra [HI UHl E] ■ This associative 
algebra T with unity consists of two generators - called standard generators - acting on a vector space V, say 
A : V ^ V and A* : V ^' V. In general, the defining relations of T depend on five scalars p, p*,7,7* and j3. In 
the following, we will focus on the reduced parameter sequence 7 = 0, 7* = and f3 = q + q^^ which exhibits 
all interesting properties that can be extended to more general parameter sequences. In this case, the so-called 
tridiagonal relations take the form 

[A,[A,[A,Al,],-i]-p[A,A*], [A*,[A*,[A*,A],],-i]=p*[A*,A], (5) 

where the g-commutator [A. 6]^ = q^'^AB — q^^'^BA has been introduced. Here, g is a deformation parameter 
assumed to be not a root of unity. Note that for p = p* = the relations ([5]) reduce to g— Serre relations, or for 
g = 1, p = p* = 16 they coincide with the Dolan-Grady relations [12] . 

The purpose of this paper is to investigate a family of TD pairs recently discovered in the context of quan- 



turn integrable systems and related algebraic structuresQ There are some reasons to do so |10| . appart from 
the applications to integrable systems that will be considered separately (TB]. Here, we study in details their 
matrix representations: eigenvalue sequences, eigenvectors and block tridiagonal stucture with respect to the dual 
basis. Using these data, we identify a family of coupled recurrence relations and a set of g— difference equations 
generalizing ([T]), ([3|), and describe their explicit solutions in the simplest cases. Our results are in agreement with 
known ones [SlITO] and provide further understanding of TD pairs. In particular, they suggest that the tridiagonal 
algebra ^ may provide a classification scheme for orthogonal symmetric functions of one argument. 

Convention: In this paper, R, C denote the field of real and complex numbers, respectively, M* — R\{0}, 
C* — C\{0}. For convenience, the deformation parameter is sometimes denoted q = exp(0) with (/> g C. According 
to a common notational convention, for a linear transformation A its conjugate transpose is denoted A*. Note 
that we do not use this convention. Here, the conjugate transpose of A is denoted A^. 

2 A family of tridiagonal pairs 

Suppose that each linear transformation A, A* is diagonalizable on V i.e V is spanned by the eigenspaces of A or 
A*. A general object called a Tridiagonal pair is defined as follows: 

Definition 2.1 Q/ Let V denote a vector space over an arbitrary field K with finite positive dimension. By a 
Tridiagonal pair (or TD pair) on V we mean an ordered pair A, A* , where A : V —> V and A* : V —^ V are linear 
transformations satisfying (i)-(iv) below. 

1. A and A* are both diagonalizable on V . 

2. There exists an ordering Vq,Vi, . . . ,Vd of the eigenspaces of A such that 

A*Vn C K-1 + K + Vn+i (0 < n < d), (6) 

where V^i ~ 0, Vd+i = 0. 

3. There exists an ordering Vq , Vi , . . . ,V^ of the eigenspaces of A* such that 

Av; c v:_ , + f; + K*+i (0 < s < <5) , (7) 

where Vl^^ = 0, V;_^^ = 0. 

4. There is no subspace W of V such that both AW C W, A*W C W, other than W ^ and W ^ V . 

We say the above TD pair is over K. 

Although a complete classification of tridiagonal pairs is not known yet, the following results have been 
obtained. Referring to the TD pair in the definition above, as shown in [^9., Lemma 4.2] one has d — 6 which is 
called the diameter. Also, for < n < d, the dimensions of the eigenspaces Vn and V* are equal and the sequence 
(dim{Vo), dim{Vi), . . . , dim{Vd)) - called the shape vector - is symmetric and unimodal i.e [T7J 

dim{Vn) — dim{Vd-n) for < n < d , dim{Vn-i) < dim{Vn) for I < n < d/2 . (8) 

Among the known examples of TD pairs, one finds a subset such that A, A* have eigenspaces of dimension one 
i.e. with a shape vector (1,1,..., 1). These are called Leonard pairs, classified in [3]. In particular, they satisfy 
(for details, see [7]) the Askey- Wilson (AW) relations ^ first introduced by Zhedanov [5]. Other examples of 



^The tridiagonal algebra with relations ||5ll is the quantum group structure behind the reflection equation associated with U 1/2 {SI2) 
I13II14I . It also plays an important role as a fundamental non-abelian symmetry in some integrable systems: it generates a hierarchy 
of mutually commuting quantities which ensure the integrability of the system. Among the examples of integrable systems which 
enjoy this symmetry, one finds for instance the boundary sine-Gordon model or the XXZ open spin chain with general boundary 
conditions. [T51 [Til ITS) . 



TD pairs are for instance the subset corresponding to p = p* = which reduce (O to the g— Serre relations of 
Ugi/2{sl2)- In this case, the shape vector associated with the TD pair A, A* is such that 

dim{Vn) < I I , 0<n<d where ( ) = ,, , ' ,, (9) 

\ n / \n I n]{d~ ny. 

denotes the binomial coefficient. It is important to mention that TD pairs also arise from irreducible finite 
dimensional representations of the Lie algebra sh and the Onsager algebra [18^. 

From now on, we focus our attention on the C/„i/2(s^2) algebra which plays a crucial role in the following 
analysis and we set IK = C The quantum Kac-Moody algebra C/^1/2 (sh) is generated by the elements {Hj, Ej, Fj}, 
j G {0, 1}. Denoting the entries of the extended Cartan matrix|2| by {dij}, they satisfy the commutation relations 

[Hi,Hj] = , [iJj, Ej] = a^jEj , [Hi, Fj] = -a^jFj , [Ei,Fj] = % — 1/2 _ -1/2 

together with the g— Serre relations 

[E,,[E„[E,,Ej]g]g-i]^0, and [F„ [i^„ [F„ Fj],],-i] = . (10) 

The sum K = Hq + Hi is the central element of the algebra. The Hopf algebra structure is ensured by the 
existence of a comultiplication A : C/ 1/2(3^2) ~* U„i/2{sl2) '^ U„i/2{sl2) and a counit £ : U„i/2{sl2) -^ C with 



(11) 

and 

£{E,)^£{F,)^£{H,)^0, £{I) = l. 

More generally, one defines the A^— coproduct A^^^ : [7^1/2(3^2) — > Uqi/2{sl2) (81 • • • (Xi Ugi/2{sl2) as 

AW = {idx---xidxA)o A(^-i) (12) 

for iV > 3 with A^^) ~ A, A^^' = id. The opposite 7V-coproduct A'(^) is similarly defined with A' = ct o A 
where the permutation map a{x (><) y) = y ^ x for all x,y €z U„i/2 (5^2) is used. 

In order to construct explicitely a family of TD pairs, we use the link recently exhibited [131 [13] between 
the tridiagonal algebra T and a class of quadratic algebra, namely the reflection equation introduced in [19l [20] . 
Based on [Ml [14], we restrict our attention to the quantum Kac-Moody algebra [7^1/2(5^2) with zero center i.e. 
K = Q. In this special case, this algebra with defining relations ([T0|) is called the quantum loop algebra of sh 
denoted [7^1/2 (£(3^2)) below. We have the following result (see also [131 114]): 

Proposition 2.2 Let {k^,k^,e±} denote nonzero scalars in C. There is an algebra homomorphism T i— > 
Uqi/2{C{sl2)) such that 

A ^ k+Eiq'^-'^ + k^Fiq^-'^ + e+q"-'^ 
A* ^ k^Eoq"°^^ + k+Foq"<'/^ + e^q""^^ (13) 



q^/^\-^-^J 




A{E.) - 


. E,®q"^^' + q-"'^'®E., 


A(F,) . 


- F,®q"^'^ + q-"''^®F,, 


A(iJ.) . 


= H,(i)I + I(i)H, 



with 



p^p*^iq'/'^+q-'/Yk+k.. (14) 



3With i,j e {0, 1}: an = 2, a^j = -2 for i ^ j. 



Proof: Replacing P^ in (O, we use the defining relations and the g— Serre relations (jlOp to simplify the l.h.s. 
of the tridiagonal relations. Straightforward calculations show that the l.h.s. of ^ reduces to the r.h.s provided 
p ^ p* is fixed in terms of the deformation parameter q, with the relation given above. ■ 

Remark 2.3 This realization is different from the one proposed by Ito and Terwilliger in JJTf in which case p = 0, 
i.e. 0) reduce to q—Serre relations. 

Remark 2.4 Using the homorphism defined above, it follows I14f that the corresponding tridiagonal relations 0) 
are invariant under the action of the coproduct ill}) as well as its generalization il2\) . 

To obtain finite dimensional irreducible representations of the standard generators (|13p . one introduces the 
evaluation homomorphism tTv ■ U^i/2{C{sl2)) i— > C/^i/2(s^2)) in the principal gradation: 

TTviEl] = VS+ , TTviEo] = VS- , 

TT^[Fi] = V^^S^ , TTviFo] == V^^S+ , 

n^q"^/^] = q^^ , Tr^q"""^^] = q-^^ , 

where the generators of Ugi/2{sl2) satisfy 

[S3,S±]=±S± and [S+,S^]^{q^-^-q-^'^)/{q'/'-q-'/'). 

As shown in |21j . every finite dimensional irreducible representation of U„i/2{C{sl2)) is a tensor product of evalu- 
ation representations. Using the A^— coproduct p2|) . it follows [T3lfT4]: 



Lemma 2.5 Let fc+ = (fc_)^,e± nonzero scalars in C. For any Vi € C* with integers i = 1, ...,N, let the pair of 
elements 1/1/q , W^ defined by 

\A/f^ = (k+v-^\-^'^S+q-'-'^ + k^Vj.q^'^S-q-'-"') ® l^^-^) + q-'^ ® W,""-'^ (15) 

with I/I/5 = e+ , Wl = e_. Let V denote a finite dimensional N~ tensor product representation of ^^ U„i/2{C{sl2)) ■ 
Assume V is an irreducible (lA^ , W\ ') — module and each of I/V5 , M^ is diagonalizable on V. Then the pair 
W^ , W^ acts on V as a tridiagonal pair. 

Proof. Conditions 1 and 4 of Definition 12. II hold by assumption. From Proposition l2.21 in the evaluation repre- 
sentation of Uqi/2{C{sl2)) the elements Wq , W^ ' satisfy the tridiagonal relations (fT^ . Using the A^— coproduct 
p2|) and the result of [21] , it follows that Wj, , Wj^ satisfy the tridiagonal relations (|13p with fM]) for general 
values of A^. Then, according to [[lOj, Theorem 10] the assertion follows. ■ 



3 Matrix representation and block tridiagonal structure 

For general values of q, irreducible (2j + 1)— finite dimensional representations of U„i/2{sl2) are well-known. In 
the following, for simplicity we will focus on the case j = 1/2 and introduce the canonical basisQ {/+,/!} for 
the i — th two-dimensional representation in the A^— tensor product representation of iS)}^U„i/ 2(312). Defining the 
endomorphism ^ : [/ 1/2(5/2) '-^ Endi^C"^), in this basis one has in terms of the Pauli matrices a±, a^: 

^+ - a+ = ( [j J ) , 5_ .- a_ = ( ; M , 253 .-. as = ( J _M . (16) 



The scalar product (., .) is defined by (/j_, /j_) = 1, (/j_, /i) = for any integer j 



Definition 3.1 Define k^ = (fc_)t — —{q^^^ — q ^^'^)e^^ /2 with 6 £ R, q — e'^ with (j) purely imaginary, and 
{a, a*} e C*. Let V = (C^)®^. Define W^°^ = e+ = cosha, wf^ = e_ = cosha*. We define the family of 
matrices of size 2^ x 2^; 



(C^)®^. 


Define W^°^ = e+ = cosh a, wf ^ = e_ = cosl 


yvr 


= (fc+a+ + A;_a_) ® J^^-^) + g-^^/^ ® W^"^"') 


w^' 


= (fc+(T+ + fc_a_) ® J(^-I) + g-'^^/2 ® WJ"^"- 



(17) 

Below, we study in details the structure of these matrices for general values of N in their eigenbasis and dual. 
Although it corresponds to the spin-j' = 1/2 case, the structure of the pair Wq , yV{ can be studied along the 
same line for higher representations leading to analogous results. 

3.1 Example A^ = 1 and the Askey- Wilson algebra 

As pointed out in 13J, for A^ = 1 the elements Wg ,W]^ ' satisfy the Askey- Wilson algebra ^ introduced in 
[5]. In particular, their matrix representations of size 2j + 1 x 2j + 1 provide an example of Leonard pairs [7] 
with shape vector (1,1,...,1). For the simplest case j = 1/2, Wq ,W} are 2x2 matrices which entries in the 
canonical basis follow from (fT7| with (fT6|) . Let ipn , n = 0,1 denote the two eigenvectors of Wq in V with 
(multiplicity free) eigenvalues A„ ' . One finds A„ ' — cosh(a + (1 — 2n)<j)/2), n = 0, 1 with 

Assuming a generic, they are independent of each other. In this new basis, the element W]; can be written as 



K' <^ 



where 



(1) cosha* sinh(Q; — (/)/2) — sinh((/)/2) (i) cosha* sinh(a + (/)/2) + sinh((?!)/2) 

sinha ' sinha 

r(i) e"(cosha + cosha*)sinh((/)/2) (i) e~"(cosha + cosha*) sinh((/)/2) 

sinh a ' sinh a 

By analogy, for generic a* € C the eigenvectors of the element W} - denoted ips with s = 0, 1 - form a 
complete basis of V. They are easily obtained substituting 

n^s, a^-a* , (j) -^ -cj) , 6 -> 6 + n (20) 

in (|18p . In this dual basis, the matrix W{ is diagonal with eigenvalues As = cosh(a* + (1 — 2s)0/2), s = 0, 1 

whereas the entries {ci^\ b, , C^^^} of W^^^ follow from {al^\ b^^\ d^^} using ^. 

3.2 Example N = 2 

Contrary to the case A^ = 1, the elements Wq , W]^ do not satisfy the defining relations ([4]) of the Askey- Wilson 
algebra [14^. Instead, they satisfy a pair of fifth order relations that can be found in [T^. Also, the eigenvalues of 

(2) (2) 

the 4x4 dimensional matrices Wq , W]; are not mutually distinct. More precisely, the shape vector associated 
with these eigenvalues A„ — cosh(a + (2 — 2n)<j)/2) for n = 0, 1, 2 takes the form {dim{Vo), dim{Vi), dim{V2)) — 
(1, 2, 1). For generic values of the parameters a, (j), 6 the corresponding eigenvectors denoted %() ,L form a complete 
basis. In the canonical basis (/|. ® /]_) of Ugi/2{sl2) <Xi Uqi/2{sl2), they take the form 

<i] = (e"+*/^+*^/U/-)®(e"+''V| + ./l), 

<i = (e-"+*/^+^VU/')«(e'"+'VU/-). (21) 



where the index i = 1,2 which characterizes the multiphcity of the eigenvalues has been introduced. With respect 

(2) (2) (2) (2) (2) (2) 

to the basis V'ofiijV'ifiii V'miJ V'2fii the matrices representing Wq ,W{ read, respectively, 



Wf =diag(A(^\Af\Af\Af), W_ 



(2) 



/ "olii] 


^(2) 

4[ii] 


^(2) 
4[12] 





7(2) 


^(2) 


"'1[12] 


^[11] 


7(2) 
f^0[21] 


"'lf211 


"'lf22l 


^(2) 
^[21] 


V 


7(2) 
"l[ll] 


7(2) 
^1(12] 


"'2[11] 



\ 



where the entries are 



a 



(2) 
'Oflll 



CLq sinha — sinh(0/2) 



sinh(a + (/)/2) 



a' 



(2) 



l(l 



ai^ sinh(a - 0) - sinh((/)/2) 
sinh(a - 0/2) ' 



(2) 



a' 



(2) 

'1[21] 



' sinh ( 



-(1) 



""2 [11] — 

f^O[ll] — 

7(2) _ 

^1(11] - 

^(2) 

'^2[11] — 

r>(2) 

^l[lll ~ 



sinh(a + </>/2) ^ ' 
a^^^ sinha + sinh(0/2) 



a: 



"'I [22] ~ 



e "sinh0 7(1) 
~sinh(a-0/2) " ' 

a[,^^ sinh(a + <j)) + smh{(j)/2) 
sinh(a + 0/2) ' 



sinh(a — 0/2) 
^0 ' 



^^^1- sinh(. + 0/2) (^° +--^(" + '^/^))' 



e"-^/^sinh(0/2) ^^(i) 
sinh(a - 0/2) 



g-a+0/2 sinh(0/2) 

sinh(a — 0/2) 
-0/2 sinh(a- 0/2) ^(1) 
sinh(a + 0/2) 

.(2) 



a\"' +cosh(a-0/2)), 

{a[^^ +cos: 



l(2) 



f^l[12] — 



, sinh(a + 0/2) ^(1) 
sinh(Q; ~ 0/2) 



^0 ' 



h(a-0/2)). 



-(2) _„4'/2pW 



-2[21] 



•>(2) 



-1[12] 



-"-^/^sinh(0/2),^(i) 



(22) 



sinh(Q; + 0/2) 



{alf' + cosh(a + 0/2)) . 



;(i) 



The set of eigenvectors ^^tL that diagonalize yV{ with s — 0,2 {k — 1) and s = 1 (fc = 1, 2) are obtained using 
(HOI) in (EH), and their eigenvalues are xi^^ = cosh(a* + (2- 2s)0/2), s = 0, 1, 2. In this dual basis the matrix W^^^ 

(2) ~ (2) l' -^ ~ (2) 

takes a block tridiagonal form similar to Wi in ([22|l . where the corresponding entries {CLg , 9^ ^ C^ } follow 

(n\ 1 (2) (n\ 

from {(4 , 0„ , C„ } using (|^ . Assuming that there are no special relations between the parameters a, 0, 6* or 
a*, 0, 61, it follows from above results that the matrices Wg , W]; satisfy the requirements of Definition [23] The 
block tridiagonal structure, which didn't appear for A^ = 1 is now clear from (j22p . In the next Section, we show 
that this property generalizes for higher values of N. 



3.3 Generalization 

To describe the family of matrices Wq , W{ in details, first we need to construct explicitely a complete basis 



;.(^ 



,{N) 



of eigenvectors ipliU] that diagonalize Wq for any N, and identify the corresponding eigenvalues. 

Definition 3.2 Let V denote the finite dimensional representation (of dimension 2^ ) on which Wq acts. 
Suppose Wq is diagonalizable on V, and let Vo,Vl,V2,... denote a standard ordering of the eigenspaces of 
Wq . For n — Q, 1,2, ..., let An denote the eigenvalue sequence of Wq . 



The structure of the eigenvectors and eigenvalues is as follows: 



Proposition 3.3 Define k^ = (fc_)' = —{q 



t = -('„l/2_„-l/2^ 



e'72. Let{a,a*} £ 



purely imaginary. Introduce the set ey — ±1, j — 1, ...,N and i G {1, ...,dim{Vn)}. The matrix Wq defined in 



be generic scalars and 

(N) 



l [J7p is diagonalized by the eigenvect 



ors 



N 



!=1 

T/ie corresponding eigenvalues, ordered by the integer n — 0,1, ..., N , are given by 

N 



Xll^) = cosh(a +{N~ 2n)(j)/2) where n = (iV - ^ ef)l2 . (24) 

fc=i 

Proof. According to previous results, the assertion is true for A^ = 1, A^ = 2. To show it for aU A^, we 
proceed by recursion. Suppose the assertion holds for A^ fixed i.e Wq VnU = ^n V'nM for a given set e' — ±1. 
Replacing A^ — > A^ + 1 in the definition pT|) and using 

CT±./=F = J± , cr±/± =0, cr3/± = ±,/± , 

we obtain 



({k+ + gl/2;^We4 + l"+4 + l ELl eL'V/2+^e) jJV+l 



(TV) 



From the definition of k± together with (|M)l and the explicit expression of n, the r.h.s. of the equation above 
reduces to 

N 

cosh (a + ( ^ 4:1 + ^15+r)0/2) (e^S^V."+4V. EL. 4-V/2-f..^A^+i + ;A^+i) ^ ^W ^ (25) 

fc=i 

provided e^ , -^ = ±1. Identifying J2k=i ^k ~ N + l — 2n, the eigenvalues take the form proposed in (|24p replacing 
A^ by Af + 1. As a consequence, the eigenvectors and eigenvalues structure for A^ also holds for A^+1. The assertion 
being true for A^ = 1, 2, it holds for all A^. ■ 

Corollary 3.4 Let {a, a*} £ C* , 6 E M. be generic scalars and (j> purely imaginary. The matrix Wq has exactly 
A^ + 1 mutually distinct eigenvalues Xk . The dimension of an eigenspace Vn is: 

dim(Vn) = ( ^ ) ■ (26) 

Proof. First, according to the definition of n in |[M|) and ej. — ±1, one has —A < X]fc=i ^L — ^- I^ yields 
< n < A^, which shows the first part of the assertion. To show ^^, we proceed by recursion. For A^ = 1, A^ = 2 
the assertion is true. Let A^ be fixed and suppose ([26]) . Notice that dim{Vn) is the total number of configurations 
with index [i] characterized by the sequence {e^ = ±1, ..., e'^ = ±1). Replacing A^ by A + 1 in ([23]) and setting 
n = (A^ + 1 — J2k=i ^fc )/2 there are two famihes of configurations associated with n: 

( (eS^l = ±1, ..., f = ±1, e3_, = Tl, ..., ^r ^1 = ±1), el^+r = +l) 
( n-1 ) configurations 
and ( (eW = ±1, ..., ^U ±1, ^^ ^ ±1, ..., ^^^ = ±1), ^^^ . -l) . (27) 

( ^ ) configurations 



It follows that Vn for iV + 1 is such that dim{Vn) = („_i) + ( „ ) = ( „ ^) i^^ agreement with ^. Then, ^ 
holds for all values of TV. ■ 

For generic scalars {a, a*, 0, 0}, the eigenvectors (j23p are linearly independent. In addition, X]n=o dim{Vn) — 
2^ = dim{V) so they provide a complete basis for V. Let us now show that yv} exhibits a block tridiagonal 
structure in the basis ^ ,-L * = 1, 2, ..., ( ). 

Proposition 3.5 Define k+ = (fc_)1^ = -{q^^'^ - g^-^/^)e**/2. Lei {a, a*} G C*, e R 6e generic scalars and (f> 
purely imaginary. Let {^„r -i }, j — 1, 2, ..., ( ^ ) as defined by \23\) the basis such that 

wfViS = ^""^^^^S ™^^ aW = cosh(a + (TV - 2n)0/2) (28) 



for n — 0, 1, ..., N . Then, the matrix W} acts on Vn as: 



(" i) (") ("-i) 

i— 1 i— 1 i— 1 

where the coefficients tt ,..,, u^t^j], C ,. , are determined recursively from il9\) . 

Proof. For N — 1,2, the assertion holds, with the entries given by (|19p and (|22l) . respectively. To show the 
assertion for all N, we proceed by recursion. Let N be fixed and suppose (|29|) is true. As noticed above, given n 
there are two families of possible configurations, given by (|?7|) . Replacing A^ by iV + 1 in ([^5)1 the corresponding 
family of eigenvectors reads: 



^Sr^ = (e"+(^-^")^/^+^^/r ^ + /^+^) V'S for je{l,...,Q}, 



n / \ n 

Having these equations, we can now consider the action of W} on the eigenvectors, using (fT7)) . Let us first 

consider the domain jG{l,...,( )}: 



(AT) 

"[i] 



("+l) 

i=l 

On the other hand, WJ; is obtained replacing A^ by A^ + 1 in (|29p . It must coincide with 

/•JV + 1\ f" + M f" + M 

-i— 1 -i— 1 i— 1 

Taking into account the two families of eigenvectors and changing appropriately the indices, (j31p becomes 

(") 

V\/^+l^//^+l) - V^ //^(A^+l)^a + (Ar-2n)0/2+ie , /,(^+l) ^ p-a-(Ar-2n)0/2+ie WA^+l 






(" 1) 

I V^ ( ( ri^N+'^^ ^-a-{N-2n+2)<jj/2+ie , p{N+l) a + (N-2n+2)(t>/2+i0\fN+l 

i=l 

, C^(W+1) , ^(W+1)n .Ar+l\ , (TV) 

i=l 

Identifying ([50)1 with (|32p . the N + 1 — th entries are determined recursively in terms of the N — th entries (see 
Appendix A) and the final result is consistent with the block tridiagonal structure. In particular, 

Then, (P5|) is satisfied for A^ ^ A^ + f and j e {1, ..., ( „ )}, provided it holds for A. We now consider the action 
of yv} ' on the eigenvectors in the domain j <={(„) + Ij ■■■, ( n^)}- We follow the same approach as above: 
after a straightforward calculation, we find similarly that (|29p for A^ ^ A^ + 1 is also satisfied, and in particular 

k,, -0 for ^e{l,...,(^_^J}, .e{(J+l,...,( ^^ )}. 

As a consequence, dH]) holds for any j G {(l, ..., (^^^). We have checked 1^ for A = 1, 2 and A' = 3 - the case 
where some entries are vanishing - so we conclude that the assertion holds for all A^. Explicit expressions for the 
entries are reported in Appendix A. ■ 

Definition 3.6 Let V denote the finite dimensional representation (of dimension 2 ) on which W} acts. 
Suppose W} is diagonalizahle on V, and let Vq ,Vi ,V2 , ■■■ denote a standard ordering of the eigenspaces of 
W} . For s — 0, 1,2, ..., let As denote the eigenvalue sequence of Wi . 



Proposition 3.7 Define fc+ = (fc_)^ = — (g^^^ — 9^-^/^)e''^/2. Let {a, a*} £ C*, 9 e R be generic scalars and 
(p purely imaginary. Lntroduce the parameters i\ — ±1, / — 1, ..., A, k £ {1, ..., ( )}. For any N, the set iPgrJ 
defined by: 



^13 = <S> (-e-^'""*-''" ^'^'^ ^U/2+^Of^ ^f._y s=iN-j: f )/2 (33) 

1=1 3=1 

form a complete basis of V . Ln this basis, the matrices Wq , W]; are such that: 

W^VIS - ^i'^ViS n^^th AW = cosh(a* + (A - 2.)0/2) , 

1=1 1 = 1 1=1 

where the coefficients Qi^nu, Ognf.T, Cguu are obtained from Appendix A with the substitution i20\) . 

Proof: According to the expressions (|17p and previous analysis, the proof is straightforward. ■ 
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Lemma 3.8 Let {a, a*} G C*, G M 6e generic scalars and (j) purely imaginary. For any N, let the matrices 
W^o'^\w[^^ be defined by (I^j. Then, W^^\ w{^^ act on V as a TD pair of diameter N. 

Proof. We show that the conditions (1-4) of Definition 12. ll are satisfied. According to Proposition 13.31 and 
Proposition 13.71 the condition 1 is satisfied. Having n G {0, l,...,Af} and the block tridiagonal structure ((29|) . 
(|34l) . V — 0„^o ^i i^ irreducible i.e. condition 4 is satisfied. Finally, Proposition 13.51 and Proposition [3l7] imply 
that the conditions (2,3) are satisfied. Then, Wq ,yV{ is a TD pair. From Corollary 13.41 and Proposition 13.71 
one finds that the diameter of the pair is d = N. ■ 



Remark 3.9 In agreement with ^9], the eigenvalue sequence (Aq , A-^ ,...,A]y ) and dual eigenvalue sequence 

(N) UN) UN 

1^1 ^■■■i^N 



Cxlr\~x[''\...rx'^^)sat^sfy 



UN) _ UN) UN) _ UN) 

'^n-2 '^n+l 1-1,1 J ^s-2 ■^s+l , -1 , -, 

T(N) — ;w='^ + '^ +^ ""^ Tw — TW='^ + '^ +^ 

/or n, s = 2, ...,7V— 1. T/ie parameters p, p* in Q) are given by p — p* — —{q — g^^)^/4. 

4 Examples of related symmetric functions 

As suggested in [TUl , given the connection between Leonard pairs|£| and the Askey- Wilson orthogonal polynomials 
(g— Racah, g— Hahn,...) on a discrete support it is interesting to investigate the family of functions associated with 
more general TD pairs, which is the purpose of this Section. As shown in [5], Askey- Wilson polynomials of general 
form with a discrete argument satisfy a three term recurrence relation and a second-order g— difference equation 
that can be derived starting from the structure of the finite dimensional representations of the Askey- Wilson 
algebra. In particular, to the four parameters of the polynomials correspond the four independent structure 
constants of the algebra. Following a similar point of view, for the TD pairs (fTS]) it is possible to deduce a coupled 
system of recurrence relations and a set of second-order q— difference equations. Below, we consider in particular 
the family ^ for a, a* purely imaginaryl so that {Ai^\ aI^^} G R. For this choice, (>V,$^V = W^^^lc-^a for 
any N. So, we define (pl^J^ = i'l[^\a^a',n~,s ■ 

4.1 Example A^ = 1 

For A^ = 1, according to previous analysis the pair Wg , yV{ is a Leonard pair. By analogy with [5], using the 
scalar product in the canonical basis /j_ we introduce the function .^„m(A) for n = 0, 1 - the index fc = 1 is 
ommited for simplicity - such that 

i^^lv^m) = u['\s):F%i\i'^) with 4llim EE 1 (35) 

for all s = 0,1. For this choice of normalization, the coefficient U{ (s) = ^a-{i-2s)a _j_ ^ ^ as follows from 
the definitions (fTS]) . Using the matrix representation p9|l for W} \ the expression {(p,L,Wl Virii) yields the 
following (two term) recurrence relations: 

Ai^^^<^i(A«) - 6i^via(A«)+a(^)^(a(A«), 

Ai^)^jfi(Ai^)) = a(VJi(A(^)) + ci^)^Ji(Ai^)) (36) 

for any s = 0, 1. For the choice ([55)l . the first equation determines uniquely !F^J,(X) = (A — (2q ) / Oq which 
form coincides with the first of the Askey- Wilson polynomials with argument A. Replacing this latter expression 



^Leonard pairs are the simplest examples of tridiagonal pairs. Their shape vector being (1, ..., 1), all eigenspaces have dimension 
one. In particular, in the basis that diagonalizes A (resp. A*), the matrix representing A* (resp. A) is irreducible tridiagonal. 
^Similar analysis can be done for a'' = a* in which case (Wg )t = W} . 
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in the second equation of (|57|) . the values of A are restricted to A = cosh(a* ± 0/2), in perfect agreement with 
the structure of the eigenvalues in the dual basis as found in Section 3.1. 

On the other hand, we may expand the expression ("^^m, Wg W'nm) using the definition (j35p . In this case, we 
now get the following (first order) g— difference equations: 

\(l)7r(l) a(l)^ - A^^^ ^l^ (l-) y(l) /x(l)x I y,(i).r-(i) n(i)^ 



A(^)^«(A«) = alVa(A«) + ci^^^^«(A«) (37) 



for any n = 0, 1. The fact that we obtain a two term recurrence relation and a first order g— difference equation 
is not surprising, as we started with spin-j = 1/2 representations of U„i/2{sl2)- Indeed, three term recurrence 
relations and second-order g— difference equations appear for spin-j > 1/2 representations of {7^1/2 (5^2) only. 

4.2 Example A^ = 2 

(2) (2) 
For N — 2, the pair Wg , Wi is the first example of TD pair that is not a Leonard pair, due to the multiplicity 

of one of the eigenvalues. We introduce the fonction ^^^Ju W for any n = 0, 1, 2 and i € {(l, ..., (^)} such that: 

i^%^^ = ui^\^)^%i^') (38) 

for s = 0, 1,2 and k £ {(l, ..., (5)}. Choosing the normalization 

.F^|^^yAf) = l for any s = 0,l,2, (39) 

and using the scalar product in the canonical basis and the explicit form of the eigenvectors (|2T|) , it follows 

C/f'(s = 0) = (e"-°'+l)2, 

C/f)(s = l) = (e"-°'+'* + l)(e"+"* +1) , C/f^(s = l) = (e"+"*+'^ + l)(e"-"*+l) , 

C/f'(s = 2) = (e"+"*+l)2. 

Then, using the matrix representation (j22p for Wj we immediately derive the following (coupled) system of 
recurrence relations - the index k is ommited for simplicity: 

^'^^(^?') - Cii^i[ii(^i'^) + di^ia(^i'') + (^n^mO^) ' (40) 

A(^)^(a(A(^)) = 6![;,^|i(A(^)) + al^j.^l^) (Ai^)) + afl^Ti^^a^P) + c(^j,^^a(A(^)) , (41) 

Ai^^^|g](^i'^) = ^![;2]^?i(Ai^^) + <]2,-^l[i(Ai^^) + <L,-^lS(Ai^)) + cl^L,^^?i(Ai^^) , (42) 

Ai^)^(f) (A(^)) . ag),^(^) (Al^)) + Cg),^lf) (Al^)) + Cg),^lf) (Al^)) , (43) 

for any s = 0,1,2. For the normalization p9p . a straightforward calculation shows that ^„rMA) are rational 
functions of A uniquely determined by the equations pp]) . (|^T|l and (|^51) for any values of A. Explicitely, together 
with ([M)l they read 

-f(2)|'i:\ _ e~"~'^/^sinh0(cosha -l-cosha*) (A - Ui[i]) 

i[i]^ ^ ^ sinh(a - 0/2) cosh a* + cosh(a + (f>/2) sinha - sinh(</«/2) (A - u) ' 

(2).r. _ e-"-'^/2gij^jj(Q,_^0/2)sinha (A - u^p])(A - ^y ) 

i[2l*- ' ^ sinh((/)/2)(sinh(a - 0/2) cosha* + cosh(a + 0/2) sinha - sinh(0/2)) (A - u) ' 

_ e-2agijji^(Q,^0/2)(cosha + cosha*) (A-U2[i]) ...^ 



_(2),-, e ^" sinh(a + 0/2) (cosh a + cosha*) (A — U2[i]) 

^[^1 ^ sinh(a - 0/2) cosh a* + cosh(a + 0/2) sinh a ~ sinh(0/2) (A - v) 
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where the expressions wiri], UiJu U2[ii and v are reported in Appendix B. Note that these expressions are invariant 
under the substitution a* -^ —a*. FinaUy, the equation (I42p fixes the values of A: replacing ^^ in (|^^ and 
using the explicit expressions of the coefficients, one gets a third order polynomial in A which roots reproduce 
exactly the eigenvalue sequence dMl) for iV = 2, as expected. 

On the other hand, it is possible to expand ((^^,y,yVg V'nM) using (|38p which leads to a set of g— difference 
equations. This being analogous to the one for higher values of TV, we refer the reader to the analysis below. 

4.3 Generalized recurrence relations and g— difference equations 

Let us now turn to general values of N. Due to the block tridiagonal structure of the matrices Wq , Wf with 
respect to their dual basis, for arbirary N it is easy to derive a system of recurrence relations and g— difference 
equations. By analogy with ((38)) . we introduce the fonctions ^^uLW for any n = 0, 1, ...,iV and i G {(l,--., („)}: 

i^M!S) = uris):Fi;i,cx^^) with 4;i,cx?^^) . i (45) 

for any s = 0, 1, ..., N and k £ {(l, ..., ( ^ )}. Here, the coefficients Ul (s) can be calculated recursively based on 
previous results: using the scalar product in the canonical basis together with the explicit form of the eigenvectors 
for s = 0, 1, ..., iV we find: 

'7V-1- 



U^is) = (e(— '+^^) + l)t/f-^)(.) for ke{l,...,{''~')}, 

U^is) = (e("+"*+(^-^)^) + l)[/(;:^-A)(.-l) for ke{{^-')+l,...,Q}. (46) 



First, let us focus on the system of recurrence relations for the rational fonctions ^^uLiX). Using the results of 
previous Sections, from {(p^J^Wi WniJ) ^^^d ommiting the index k we immediatly get 

^^f(^) - E d^:i,](A) + E<3,^HA)+ g CS)^-),,(A) (47) 

2—1 i—1 i— 1 

for A e {\s }, n,s — 0, 1, ..., N, j e {1, ..., {n)} with the coefficients given in Appendix A. Notice that similarly 
to the case N — 2, 2^ — 1 equations will determine uniquely the explicit form oi T^^^J {X) . The remaining equation, 

associated with n = A^ — 1 and j = N, determines the eigenvalues to be in the set A G {As }, s = 0, 1, ...,N 
given by ((M|) . 

To derive asystemofg— difference equations for the rational fonctions ^nuLiXs), we start from {'fg\J,y^o '^nkU 
Let us introduce the discrete g— difference operator rj g{Xs) — (?(As±i). Combining the normalization coefficients 
ul (s) given by (|46p with the "dual" entries we obtain a degenerate family of second-order g— difference operators 
such that 

dW(.)f2!(AW) = Ar)F2i(AW) with F^l ^ [^S,^S, -.^IS"/ (48) 

where 

D(~)(5) = a>W(s),7 + ¥^'^'(s),7-i + ^^^^Hs) , 



ijfw ' "^''^ ''''~^Jri^ 



Ul (s) 
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It should be pointed out that for A^ = 1, the operator D(^'(s) is a special case (as we choosed the spin— 1/2 
representation of C/ 1/2(3^2)) of the Askey- Wilson operator derived in [S]. As a consequence, for higher values of 
N the q— difference operators D*^^-' (s) can be seen as some generalizations of the Askey- Wilson operator. Although 
it is an interesting problem to explore some consequences in integrable systems and the links with a recent work 
on Bethe ansatz and 5— Sturm-Liouville problems [22], we do not pursue the analysis further and leave it for a 
separate work |16j. Note that for iV = 2 we have checked explicitly that the rational functions (|44|) satisfy the 
degenerate set of g— difference equations (|48|) . 

4.4 Orthogonality 

The system of Askey- Wilson polynomials is orthogonal with respect to a certain discrete weight function. This 
can be shown, for instance, using the completness of the two dual basis which diagonalize the elements of the 
Askey- Wilson algebra [5] . Such property can be extended to the system of rational functions !F^, .^ (A) (defined on 
a real discrete support) obeying (|T7|) . ([35]), having in mind that the two basis ([^ . ([55)1 are complete for general 
values of the parameters. To show that, we introduce the basis of dual eigenvectors <^ rJ = V" ih]\a^a*,n^s and 



V'im = f„ii]\a'^a,s^n- It is casy to check that 



r 


..o,...,("-),, 




for 


..,(->,.. 


■•(" 

V n 



where, using the explicit expression of the eigenvectors, one can easily obtain recursively the normalization 
coefficients. For instance, for n = 0, ..., iV one finds 

A/"^^^ = /V(^-l)('l_e(2a+(JV-l-2n)0)N-l 
n[i] n[i] \ ) 

<• - Aff-.'V-,,(1 -.'--->"«-'«)-' for i.((";') + l,„^.(;')}, (50) 

and similarly AQr^i — N^\J\{n,i.a,<ti}^{s.k.-a'.-<ji} for s = 0, 1, ...,N. Then, the completness of the two eigenbasis 
together with the definition of the normalization coefficients yield 

s=0 fc=l 

Replacing (|45p in the above equation, and using ('4'nii] ,'-Ps\kV ~ ('^slkl^^mlV ^^ ^^'^ ^-'^^^^ *^^ system of rational 
functions J-"^.y(As) of the discrete argument A^ are orthogonal on the A^ + 1 points on the interval of the real 
axis < s < N. The condition of orthogonality reads 

E E -r(^)^S(^^)^iM(^«) - i^SV'KrnS., where w^^\s) ^ M^S {U^ is)f 

s=0 fe=l 

is the discrete weight function that ensures the orthogonality of the sytem. Given N, its explicit expression is 
derived recursively using ([50]) . ([^5]) . In particular, it is an exercise to check the orthogonality condition for A^ = 1 
(in which case our construction reduces to a special case (spin— 1/2) of [5J) and N — 2. 
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Appendix A: CoefRcients of the matrices An , Bn , Cn 

We can organize the set of eigenvectors by eigenspaces Vn and set ^„ = [V'nfij i '0„[2ii ■■■' '/'l[r«)i] ■ Then, the 
relations ((28)) . ((29|) can be written in the form 

^W^W ^ sr)<\+AW*W+CWv,W , (51) 

where An ,Bn ,Cn are submatrices with entries {Xn )iji] = (2^ )[y] which have been determined recur- 
sively in Section 3, with "initial" conditions (|19p . The size of each submatrix is determined by the multiplicity of 
the eigenvalue. We obtain: 



Matrix aI^^ 



/^-lU . nW - <r' ^^"^(" + (JV - 2 - 2n)cb/2) sinh(0/2) 



For i 6= il P-^U • a^ ■• - 

jfcli,...,^„ ;|. a„[^,,]- sinh(a + (TV - 1 - 2n)</)/2) 

af";'^ .«-l^ ^ .«-l^, sinh(a + (TV + 2 - 2n)(b/2) + sinh(0/2) 
J^or ze{(„ j + l,...,(„j|. a„[^^, - sinh(a + (TV + 1 - 2n)0/2) 

17 -^/i ^^-lU ■^S(N-i\^, ^wu nW e-"-(^-^")^/^sinh(0) ^(JV-i) 

For z £{!,...,(„ )},,e{(„ )+!,...,(„)}: <[.jl = ' 3i,h(a + (iV - 1 - 2n)0/2) ^-^faH"-)! 

For *e 1(^-1) + 1 (^)l 7611 (^-1)1 ■ a(^) - e"+^^-^"^^/^sinh(0) (^,,) 

For ziGll f^-in i^7 • ^(iV) _ sinh(a + (JV - 2 - 2n)0/2) ^n^i) 

For Z7eif^-M + 1 f^n z^7- ^(JV) _ sinh(a + (iV + 2 - 2n)0/2) (Ar_i) 

J^or z,jeiU j + i.-'UJi' *7^J- "»M-sinh(a+(7V + l-2n)<^/2) "-i[-(^-^):^-(^-^)l ■ 

• Matrix B^^'' 

For .e{l,...,C-i)}, ,e{l,...,(r^)}: 6i;;].j=e^/^62,:j'\ 

For ze{l,...,(^-i)}, ,e{(r^) + l,. ..,(::)}: ^i^;;]] = , 

For *e{&l^) + l,...,&0}- ^•e{(r^)+l,...,(^)}: 



r (iV) _^/2sinh(a + (iV + 1 - 2n)0/2) , (W-i) 

^"M=^ c^.h(a + (TV - 1 - 2n)<^/2) ^"-^I-^-/) ^-(^"^)l 



sm. 



For *e{e+7) + l,...,C^+0}: 

^"I» -(".-)] = sinh(.+ (A.-l-2nW2) (<t;b) ^-(:.7)] + ^°^^(" + ^^ " ^ " ^^^^Z^)) ' 

For *e {(^+7) + l,...,(^+OLj-e {!,..., (r^)}, ^^J + (^+7): 

7W _ e°+(^-i-^")^/^sinh(0/2) (a,_i) 
"I'^l sinh(a + {N-1- 2n)0/2) »[*~("+7) j) ' 

• Matrix ck^'> 

For *G{C-r/) + l,-.(^i)}, je{l,...,Gr^)}: ci"?]] = , 
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For ie\l r-^U 7 6 11 P-^U- c^^^ - r-^/^ '^"^^" + (A^ - 1 ~ 2n)0/2) (jy^i) 



"fel ' sinh(a + (iV + 1 - 2n)0/2) "^1 



For ze{l,...,(^_7)}: 



7 

For ze{l,...,(^-,i)},je{(ri) + !,...,(::)}, ^7^J-(r^): 
W e-°-(^+i-^")^/^sinh(0/2) (^_i) 

"fol sinh(a + (iV + 1 - 2n)0/2) n-ifo-(^-i)] " 



Appendix B: Zeroes and poles of the rational functions for N = 2 

Ui[i] = — cosha , 

_ sinh(a + 2a* + 0/2) + sinh(a - 2a* + 0/2) + smh(2a + a* + 30/2) + sinh(2a - a* + 30/2) 
^ 2(sinh(a - a* - 0/2) + smh(a + a* - 0/2) + smh(2a + 0/2) - 3sinh(0/2)) 

sinh(a + 0/2) - 3 sinh(a - 0/2) + sinh(a - 30/2) - 3 sinh(a + 30/2) + 3 sinh(a* + 0/2) - 3 sinh(a* - 0/2) 
2(sinh(a - a* - 0/2) + sinh(a + a* - 0/2) + sinli(2a + 0/2) - 3sinh(0/2)) ' 

(±) [/±2sinh(0/2)\/7 

i[2l 4(cosh(0/2)-cosh(2a + 0/2)) ' 

_ sinh(a - a* - 0/2) + sinh(a + a* - 0/2) - sinh(0/2) - sinh(30/2) 
"^[^1 " 2sinh(a + 0/2) ' 

where 

[/ = 4 cosh(a + 0/2) + 2 cosh(a* + 0/2) + 2 cosh(a* - 0/2) - 2 cosh(a + 30/2) - 2 cosh(a - 0/2) 

- cosh(a* + 2a~ 0/2) - cosh(a* - 2a + 0/2) - cosh(a* +2a + 30/2) - cosh(a* - 2a - 30/2) , 

V = (cosh(2a) + cosh(2a + 0) - 2) x 

X (8 + 4cosh(a + a*) + 4cosh(a - a*) + 4cosh(a + a* + 0) + 4cosh(a - a* + 0) + 4cosh(0) - 2cosh(2a) 
+2 cosh(2a + 0) + cosh(2a + 2a*) + cosh(2a - 2a*) + cosh(2a + 2a* + 0) + cosh(2a - 2a* + 0)) . 
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